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Abstract 

We consider the problem of maximizing the revenue raised from tohs set on the arcs of 
a transportation network, under the constraint that users are assigned to toll-compatible 
shortest paths. We first prove that this problem is strongly NP-hard. We then provide a 
polynomial time algorithm with a worst-case precision guarantee of ^ log2 niT + l, where 
rriT denotes the number of toll arcs. Finally we show that the approximation is tight 
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with respect to a natural relaxation by constructing a family of instances for which the 
relaxation gap is reached. 

Keywords: network pricing, approximation algorithms, Stackelberg games, combi- 
natorial optimization, NP-hard problems. 

1 Introduction 

This paper focuses on a class of bilevel problems that arise naturally when tariffs, tolls, or 
devious taxes are to be determined over a network. This class of problems encompasses 
several important optimization problems encountered in the transportation, telecommu- 
nication, and airline industries. Our aim is twofold: first, we show that the problem is 
NP-hard; then we present a polynomial time algorithm with a tight worst-case guarantee 
of performance. 

Bilevel programming is a modelling framework for situations where one player (the 
"leader") integrates within its optimization schedule the reaction of a second player 
(the "follower") to its own course of action. These problems are closely related to 
static Stackelberg games and mathematical programs with equilibrium constraints (or 
MPECs, see Luo, Pang and Ralph ^2]), in which the lower level solution characterizes 
the equilibrium state of a physical or social system. Bilevel programs allow the modelling 
of a variety of situations that occur in operations research, economics, finance, etc. For 
instance, one may consider the maximization of social welfare, taking into account the 
selfish behavior of consumers. It is well-known that the taxation of resources and services 
at marginal cost (Pigovian taxes, see Pigou ^Hl) maximizes global welfare. However, 
when some resources fall outside the control of the leader, the social optimum might not 
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be reachable, yielding a "second-best" problem of true Stackelberg nature (see Verlioef, 
Nijkamp and Rietveld JSj, Hearn and Ramana and Larsson and Patriksson [S] for 
traffic examples). In constrast with these studies, we adopt the point of view of a firm 
involved in the management of the network but oblivious to social welfare; the firm's 
only goal is to maximize its own revenue. 

Bilevel programs are generally nonconvex and nondifferentiable, i.e., to all practi- 
cal extent, intractable. In particular, it has been shown by Jeroslow [Jj that linear 
bilevel programming is NP-hard. This result has been refined by Vicente, Savard and 
Jiidice |2(Jj , who proved that obtaining a mere certificate of local optimality is strongly 
NP-hard. Actually, Audet et al. ^ unveiled the close relationship between bilevel pro- 
gramming and integer programming. This "intractability" has prompted the develop- 
ment of heuristics that are adapted to the specific nature of the instance under con- 
sideration, together with their worst-case analysis. Such analysis was first performed 
for a network design problem with user-optimized fiows by Marcotte JHl) who proved 
worst-case bounds for convex optimization based heuristics. More recently, worst-case 
analysis of Stackelberg problems has been applied to job scheduling and to network 
design by Roughgarden |161 117j . to network routing by Korilis, Lazar and Orda [S] and 
to pricing of computer networks by Cocchi et al. All these works focus on "soft" 
Stackelberg games, where the objectives of both players are non-confiicting, and where 
heuristics are expected to perform well in practice, although their worst-case behavior 
may turn out to be bad. 

In this paper, we analyze an approximation algorithm for the toll optimization 
problem (MaxToll in the sequel) formulated and analyzed by Labbe, Marcotte and 
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Savard jlUj . In this game, which is almost zero-sum, a leader sets tolls on a subset 
of arcs of a transportation network, while network users travel on shortest paths with 
respect to the cost structure induced by the tolls. Labbe et al. ^] proved that the 
Hamiltonian path problem can be reduced to a version of MaxToll involving negative 
arc costs and positive lower bounds on tolls^. In this paper, we improve this result by 
showing that MaxToll, without lower bound constraints on tolls, is strongly NP-hard. 
Next, in the single-commodity case, we provide a polynomial time algorithm with a per- 
formance guarantee of ^ log2 rriT + 1, where mT denotes the number of toll arcs in the 
network. We then use this result as well as specially constructed instances to prove the 
tightness of our analysis, as well as the optimality of the approximation factor obtained 
with respect to a natural upper bound. 

The rest of the paper is organized as follows. In Section 121 we state the problem and 
prove that it is NP-hard. In Section |21 we introduce an approximation algorithm whose 
performance is analyzed in Section^ 

2 The model and its complexity 
2.1 The model 

The generic hilevel toll problem can be expressed as 

max Tx 

T 

It was also shown recently by Marcotte et al.^l] that the TSP is a special case of MaxToll. 
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where x is the partial solution of the parametric linear program 



mm 



(ci + T)x + Cay 



s.t. 



Aix + A2y = h 



x,y >0 



In the above, T represents a toll vector, x the vector of toll commodities and y the 
vector of toll-free commodities. 

We shall consider a combinatorial version of this problem. Let G = {V, A) be a 
directed multigraph with two distinguished vertices: the origin s (zV and the destination 
t € V. The arc set A is partitioned into subsets At and Ajj of toll and toll-free arcs, of 
respective cardinalities rriT and mu. Arcs are assigned fixed costs c : At N™^ and 
d : Ajj — N*"*^ (in the sequel, N = {0, 1, . . .}). Once tolls are added to the fixed costs 
of At, we obtain a toll network N't = {G, c + T, d, s, t). Denoting by SVIMt] the set 
of shortest paths from s to t, we can then formulate MaxToll as the combinatorial 
mathematical program (see also Figure ^): 



This is a single-commodity instance of the toll setting problem analyzed by Labbe et 

ai. unj. 

In this framework, the leader must strike the right balance between low toll levels, 
which generate low revenue, and high levels, which could also result in low revenue, 
as the follower would select a path with few toll arcs, or even none. An instance of 
MaxToll is iUustrated in Figure |2 

Several remarks are in order. Firstly, to avoid a trivial situation, we posit the 



max > 

T>0 ^—^ 




(1) 



5 



INSTANCE: 


- a directed graph G = {V, At U Au) 




- fixed cost vectors c : At ^ N"^^ and d : Au ^ N™^ 




- distinguished vertices s,t &V such that there 




exists a simple path from s to t in [V, Au) 


SOLUTION: 


- a nonnegative toll vector T on At 




- a simple s — t path P of minimal length w.r.t. 




the cost structure (c + T, d) 


MEASURE: 


- maximize EeeA^nP T{e) 



Figure 1: MaxToll 



existence of at least one toll-free path from s to t. Secondly, our formulation implies 
that, given ties at the lower level (the user's level), the leader chooses among the toll- 
compatible shortest paths the one travelled by the follower. Note that a risk- averse 
leader could always force the use of the most profitable path by subtracting a small 
amount from every toll on that path, thus yielding a revenue as close as desired to 
the revenue generated by the "cooperative" solution. Thirdly, once a path P has been 
selected by the leader, toll arcs outside P become irrelevant. In practice, the removal 
of these arcs can be achieved by setting tolls to an arbitrarily large value on toll arcs 
outside P. We denote by Mt{P) the network where these arcs have been removed. 
Finally, our central results hold also in a version of MaxToll where T is unconstrained; 
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Figure 2: This network contains two toll arcs (represented by dashed 
arcs). Fixed costs are given by numbers close to each arc. The optimal 
path is (s, V2, v^, v^, t) with a revenue of 4 when Ti = 2 and T2 = 2. 

in this case, negative tolls can be interpreted as subsidies. Actually, Labbe, Marcotte 
and Savard HjO have constructed instances where the optimal solution involves negative 
tolls. Nevertheless, throughout most of the paper we focus on nonnegative tolls because 
(i) this case is interesting in its own sake, (ii) intermediate results are easier to interpret 
when tolls are thought to be nonnegative. 

A natural upper bound on the leader's revenue has been derived by Labbe et al. ^0] 
using duality arguments from linear programming theory. It also follows from Theorem|l] 
of Section l3.ll Let C^o be the length of a shortest toll-free path and let C{P) be the 
length of a given path P with T{e) = for all toll arcs e. 

Theorem 1 Let P be a path. Then the optimal revenue associated with P is bounded 
above by 

B{P) = Coo- C{P). (2) 
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Since does not depend on P, it follows that the largest upper bound corresponds to 
the path with smallest value of C{P); that is, P is a shortest path when tolls are set to 
0. We denote the length of such a path by £0 and by 

LP = Coo- Co (3) 

the value of a path-independent upper bound. This bound is simply the difference 
between the costs of shortest paths corresponding to infinite and null tolls, respectively. 
Note that, if the set of toll arcs is a singleton, the upper bound can always be achieved. 

2.2 NP-hardness of MaxToll 

The purpose of this section is to show that MaxToll is strongly NP-hard. We also 
prove that a version of MaxToll where the toll vector is unconstrained shares this 
property, thus settling a conjecture of Labbe et al. ^0] about the complexity status of 
the generic toll setting problem. 

Theorem 2 MaxToll is strongly NP-hard. 

Proof: Let C denote the sum of all fixed costs. It is not difficult to show that there 
exists an optimal toll vector T that is integer- valued and less than C + 1; in particular, 
optimal solutions are of polynomial size. 

Now, consider a reduction from 3-SAT to MaxToll (see Let xi, . . . , Xn he n 
Boolean variables and 

m 

F= V/i2 VZi3) (4) 

i=l 

be a 3-CNF formula consisting of m clauses with literals (variables or their negations) 
lij. For each clause, we construct a sub-network comprising one toll arc for each literal 
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Figure 3: Sub-network for clause {In V li2 V k^). 

as shown in Figure 01 

The idea is the following: if the optimal path goes through toll arc Tij, then the 
corresponding literal kj is true (note: if kj = Xk, then Xk =false). The sub-networks 
are connected by two arcs, a toll-free arc of cost 2 and a toll arc of cost 0, as shown in 
Figure |31 

If F is satisfiable, we want the optimal path to go through a single toll arc per 
sub-network (i.e., one true literal per clause) and simultaneously want to make sure 
that the corresponding assignment of variables is consistent; i.e., paths that include a 
variable and its negation must be ruled out. For that purpose, we assign to every pair of 
literals corresponding to a variable and its negation an inter-clause toll-free arc between 
the corresponding toll arcs (see Figure 0}. As we will see, this implies that inconsistent 
paths, involving a variable and its negation, are suboptimal. 

Since the length of a shortest toll-free path is m + 2{m — 1) = 3m — 2 and that of a 
shortest path with zero tolls is 0, 3m — 2 is an upper bound on the revenue. We claim 
that F is satisfiable if and only if the optimal revenue is equal to that bound. 
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Figure 4: Network for the formula (xi V X2 V X3) A (x2 V X3 V X4) A (xi V 0:3 V x^). 
Inter-clause arcs are bold. Path through T12, T22, T^2 is optimal {x2 = =true). 
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Assume that the optimal revenue is equal to 3m — 2. Obviously, the length of the 
optimal path when tolls are set to must be 0, otherwise the upper bound cannot be 
reached. To achieve this, the optimal path has to go through one toll arc per sub-network 
(it cannot use inter-clause arcs) and tolls have to be set to 1 on selected literals, C -|- 1 on 
other literals and 2 on tolls T^, V k. We claim that the optimal path does not include a 
variable and its negation. Indeed, if that were the case, the inter-clause arc joining the 
corresponding toll arcs would impose a constraint on the tolls between its endpoints. 
In particular, the toll immediately following the initial vertex of this inter-clause arc 
would have to be set at most to 1, instead of 2. This yields a contradiction. Therefore, 
the optimal path must correspond to a consistent assignment, and F is satisfiable (note: 
if a variable and its negation do not appear on the optimal path, this variable can be 
set to any value). 

Conversely if F is satisfiable, at least one literal per clause is true in a satisfying 
assignment. Consider the path going through the toll arcs corresponding to these literals. 
Since the assignment is consistent, the path does not simultaneously include a variable 
and its negation, and no inter-clause arc limits the revenue. Thus, the upper bound of 
3m — 2 is reached on this path. 

Finally, note that the number of arcs in the reduction is less than 

10m + 2(m + {3mf 
and that all constants are polynomially bounded in m. 

□ 
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It is not difficult to prove that the same NP-hardness reduction works when negative 
tolls are allowed. 

Theorem 3 MaxToll is still strongly NP-hard when negative tolls are allowed. 

Proof: We use the same reduction as in the nonnegative case. The proof rests on 
two results proved in ^0]. First, the upper bound is valid when T is unrestricted. 
Second, there exists optimal solutions of polynomial size. The latter result follows from 
a polyhedral characterization of the feasible set. 

From the first result, we know that the 3m — 2 upper bound on the revenue is 
unchanged. On the other hand, if F is satisfiable, the feasible solution considered in the 
nonnegative case still yields a 3m — 2 revenue. We only have to make sure that negative 
tolls cannot produce a 3m — 2 revenue when F is not satisfiable. Again, to reach the 
upper bound, one has to use a path of length when tolls are set to 0. Consequently the 
optimal path comprises exactly one literal per clause. Now the toll-free arcs of length 1 
and 2 limit the Tj^'s on the path to 1 and the T^'s to 2, so negative tolls are useless in 
this case. Indeed, inconsistent paths will see their revenue limited by inter-clause arcs 
without any possibility to make up for the loss incurred from negative tolls. 

□ 

3 Approximation algorithm 

In this section, we devise a polynomial-time approximation algorithm for MaxToll. 
Such an algorithm is guaranteed to compute a feasible solution with objective at least 
OPT/a, where OPT is the optimal revenue and q, which depends on the number 
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rriT of toll arcs, denotes the approximation factor. For a survey of recent results on 
approximation algorithms, the reader is referred to Hochbaum jH], Vazirani |18j and 
Ausiello et al. 

3.1 Preliminaries: characterizing consistent tolls 

The leader is only interested in paths that have the potential of generating positive 
revenue. This remark warrants the following definitions. Recall that J\[t{P) is the 
network Mr in which toll arcs outside the path P have been removed. 

Definition 1 Let denote the number of toll arcs in a path P from s to t . We say 
that P is valid if > 1 and P is a shortest path with respect to a null toll vector T , 
I.e., P G SV[J\fo{P)]. 

It is clear that non-valid paths cannot generate revenue. Any valid path P can be 
expressed as a sequence 

P = {vo,l, Ti, Vi,2, T2, ■■■,T^P, f„^^;,mf+l) (5) 

where Tj is the i-th toll arc of P (in the order of traversal) and fj^j+i is the toll-free 
subpath of P from the terminal vertex TERM(rj) of Tj to the initial vertex iNiT(Tj+i) 
of Tj+i. According to this notation, uo,i starts in s and ^^m^ m^+i ends in t. Since 
P € SV[Afo{P)], Vi^i+i is a shortest toll-free path from TERM(rj) to iNiT(ri+i). We 
extend this notation to Vij, a shortest toll-free path from TERM(rj) to iNiT(rj), with 
length Uij. For k < I, let Ck,i be the length of P (i.e. the sum of costs) from TERM(Tfc) 
to iNiT(Ti) with tolls set to 0, and Tk^i be the sum of tolls between TERM(rfc) and iNiT(r/) 
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on P, 

l-l l-l l-l 

Ck,i = ^Ui^i+i+ ^ c{Ti) Tk,i= ^ T{Ti), (6) 

i=k i=k+l i=k+l 

with the convention that Yl\=k Xi = ii I < k. 

Definition 2 Let P be a valid path. The toll vector T is consistent with P if P ^ 
SVlNriP)]; that is, P remains a shortest path when tolls outside P are removed and 
tolls on P are set according to the vector T . 

The following result, which characterizes consistent tolls, is the starting point of our 
algorithm. 

Theorem 4 Let P he a valid path. Then, the toll vector T is consistent with P if and 
only if 

A,j + Tij < Uij V < i < j < + 1. (7) 

Before providing the proof of this theorem, we give an example of its application. Con- 
sider again the network of Figure El Take the path P = (s, V2, V3, V4, t). Using the 
renumbering system introduced above, we have that ri = (s, ^2), T2 = (f4, t), vq^i is 
empty (it actually consists in the node s by itself), vi^2 = (^2, ^^3, ^4) and ""2,3 is also 
empty. Other toll-free subpaths that play a role in Theorem 0] are 1)0,2 = (s, vi, V4), 
vi,z = {v2, V3, V5, t), and ^0,3 = {s, vq, t). The corresponding lengths are ^0,1 = ^1,2 = 
^2,3 = 0, Z^o,2 = 11) ^1,3 = 7, and Z//o,3 = 4. Thus, the inequalities in Theorem |1] are 

Tin) < 11, T{t2) < 7, Tin) + T{t2) < 4. 

From this we conclude that one possible set of optimal taxes is T{ti) = T[t2) = 2 as 
claimed in Figure [2 Figure [5] provides an equivalent representation of this example 
according to Theorem^ 
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Figure 5: Equivalent representation of the example of Figure |21 AH 
arcs have been replaced by shortest paths between toll arcs and ori- 
gin/destination. 

Proof of Theorem^ Obvious from the very definition of consistent tolls. <^^) 
We need to check that the length of P remains smaller than or equal to the length of 
any other path when tolls satisfy Q. This looks rather obvious. However, one must be 
careful about paths that borrow a subset of the toll arcs of P, possibly in a different 
sequence (note that toll-free paths are taken care of by setting i = and j = mf^ + 1 
above). Assume, by contradiction, that such a path P is strictly shorter than P in 
MriP), and that conditions Q are satisfied. We note 

P = {vo,i, h, Vl,2, . . . ,f^P, ^P+i) (8) 

with corresponding Uij, Ck^i and Tk^i for all i, j, k, I with k < I. For all i, fj = Ts(^i) for 
some injective function 6 (toll arcs outside P are irrelevant). To derive a contradiction, 
we will construct a path that is shorter than P by modifying the "backward" toll-free 
subpaths of P. This improved path will happen to be P. See FigurelHlfor an illustration. 
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V0,1 V2,3 




Figure 6: The thick path is P while the straight thin path is P. An 
example of the "backward" toll-free subpaths referred to in the proof is 
51,2 T2 ^ V2,3- This gets replaced by the section of P between Ts(ji) 
and Ts(j^), which is clearly shorter. 

Let ji = 1 and jk = min{j > j^-i : 5{j) > as long as such jk exists, the last 

one being denoted jx- We further define (5(0) = jo = and 6{m^ + 1) = JK+i = f^j- + 1- 
We claim that 

^Hjk-i),sUk) + '^{ife-i).'5(ifc) — ^jk~i,jk +'^fc-ijfe' 

for all k. This implies that we can replace the subpath of P between TERM(fjj,^) 
and iNiT(fj,,.), by the subpath of P between TERM(T5(j^._j)) and miT^Tg^j^^) without 
increasing its length. But after doing this for all k, we obtain P which is a contradiction. 
We divide the claim above in two cases. Firstly, if jk — ^ = jk-i (i-e. P "goes in the 
direction of P", e.g. -uq,! in Figure inj, then 

where we have used ((TJ (note that in this case, there are no toll arcs on P between 
TERM(fjj._^) and iNiT(fjj.)). Otherwise, jfe — 1 7^ jk-i (i.e. P "goes in the direction 
opposite to P", e.g. vi^2 — > T2 — J- ^2,3 in Figure IHI). By definition of jk, we have 
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'^(jfc~l) ^ <^(jfc-i) for all 1 < A; < K+l. For examph 



le, in Figureini we have ts^j^-i) = T2 



which clearly comes before ^^(ja) = "^3 -P- From this and ((TJ, we get 



□ 



Remark. There always exists an optimal solution with tolls less than C + 1 (see 
Theorem 121 for a definition of C). Indeed, Uij < C + 1, Vi, j implies that optimal tolls 
on P have to be lower than C + 1. Therefore, fixing tolls to C + 1 outside P generates 
no additional active constraints on the tolls of P. 

3.2 High-level description of the algorithm 

ExploreDescendants is an approximation algorithm motivated by the previous char- 
acterization of consistent tolls. Initially, the algorithm computes the optimal revenue 
associated with a shortest path in SV[Mo], i.e, a path with the largest upper bound 
B{P)'^. If revenue is smaller than the upper bound LP, Theorem implies that there 
exists a toll-free subpath Vk,i k < I whose short length forces some tolls in P to be small. 
To relax this constraint, it makes sense to skip the subpath of P between TERM(Tfc) and 
iNiT(r/) and to replace it by Vk,i- This yields a new path whose length will not be much 
larger than the length of P, and for which some constraints ^ have been removed. 

This path is a natural candidate for improved revenue. By repeating this process, we 
^As shown in the next section, this can be achieved in polynomial time. 
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will show that Algorithm ExploreDescendants can uncover, in polynomial time, a 
path with good approximation properties. 

Algorithm ExploreDescendants is streamlined in Figure [Tj It comprises two 
subroutines, MaxRev and TollPartition. MaxRev computes the largest revenue 
compatible with the shortest path status of P. Starting from P, TollPartition gen- 
erates two descendants of path P. The algorithm is initialized with Pq in A/q, a shortest 
path of length Cq. 
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Algorithm ExploreDescendants 
Input: a path P 

Output: a path P, toUs T and objective value V 

• Compute maximum revenue Vp achievable on P and corresponding toll vector 
Tp: {Vp,Tp,{{i'{k),f{k))}'^i^ :=MaxRev(P) 

• If Vp < B{P) then 

— Derive new paths from P : 

(Pi, P2) := TollPartition (p, /(A;))}^i) 

— For 2 = 1,2 : {Vi, Ti, Pi) := EXPLOREDESCENDANTS(Pi) 

— Return V := max{yp, Vi, V2} and corresponding toll vector T and path 
P 

• Else return (F, T, P) := (Vp, Tp, P) 

Figure 7: Algorithm ExploreDescendants 
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3.3 Maximizing path-compatible revenue 



Let P be a valid path in Af denoted by 



P = {^0,1, Tl, 1^1,2, ■ ■ ■ 




(9) 



with corresponding values olUij, C^^i and 7^ ^ for all i, j, k, I, with k < I. The optimal 
tolls compatible with P's shortest path status can be obtained from a simple greedy 
algorithm (see Figure |H1) : consider each toll arc in order of traversal, and fix its value 
to the largest value allowed by Theorem |3 taking into account the tolls previously set. 
This leads to the recursion 



The validity of the above formula rests on Theorem [S] where we construct a sequence of 
toll-free subpaths such that (i) every toll of P is bounded from above by at least one 
path in the sequence (condition below), (ii) the sum of the tolls defined by (|1U() is 
equal to the sum of the bounds imposed by the paths of the sequence (condition (|14|) 
below). 

Labbe et al. ^0] give another polynomial time algorithm for this task but it does 
not provide the information we need regarding the active toll-free subpaths. This in- 
formation will be instrumental in generating new paths from P and in obtaining an 
approximation guarantee. 

Theorem 5 Let Vij, < i < j < m^ + l, be the toll-free subpaths defined in Section Wl\ 
and let tji, 1 < k < be as in !ll(J\) . Then, there exists a sequence of paths 




(10) 



(11) 
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Algorithm MaxRev 
Input: a path P 

Output: tolls Tp, objective value Vp, and sequence of indices {(^'(^), 

• Denote P = {vo,i, n, Vi^2, T2, . . . , t^p, v^p^^p^^), and set k := 1 

• While k < + 1, do 

— Compute 

{k-l 
1=1+1 

— Let {i'{k), j'{k)) be the pair of indices for which the minimum above is 
attained, breaking ties by selecting the largest j'{k) and corresponding 
smallest i'{k) 

- Set {i'{l),f{l)) := {t'{k),j'{k)) for all A; < / < f{k) 

- Set k := j'{k) 

• Let Tp be +00 outside P and Vp := J2k=i'^k 

• Return Tp, Vp and {(^'(A;), 

Figure 8: Algorithm MaxRev 

with i(l) = and j{q) = ?ny + 1 such that for all h (see Figure\^ 

i{h + l) < j{h) <i{h + 2) < j{h + l) (12) 
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and for all k 

\{h : iih) + 1 <k <j{h) -1}\>1 (13) 
with equality iftk 7^ 0. This implies that for all h', h" with h' < h" 

j(h")-l h" 

tk=Yl [^mjw - • (14) 

k=i{h')+l h=h' 

Corollary 1 The toll assignment defined by lll(J\) is optimal for P. 

Proof of Corollary^- For every toll assignment T' consistent with P 

p 

rrirp q 

Y,T'{Tk) < by and CI 

k=l h=l 

P 

= J2^k, bydll. 

k=l 

□ 

Proof of Theorem\^ This proof is rather cumbersome, so we start by giving an intuitive 
description. By Theorem |1J we know that tolls are constrained by toll-free subpaths. In 
order to obtain the corollary above - which is the ultimate purpose of this theorem - 
we need to find a subset of constraints such that each toll arc is covered by at least one 
constraint and such that the sum of the slacks of the constraints is the same as the sum of 




Figure 9: A section of a path P, showing toll arcs with k = i{h + 
l),j{h),i{h + 2),j{h + l). 
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the tfc's. Natural candidates for saturated constraints are derived from (jlUp : when we fix 
tk, at least one constraint becomes saturated. These constraints will be denoted by the 
sequence of indices {{i' (k) , j' (k))} ^^^^ below. The problem with them is that they might 
be redundant. For example, in setting ti we might saturate the path ^0,3 and then in 
setting ^3 we might saturate the path vq^^ which "supersedes" ^0,3. This situation makes 
it impossible to obtain equality in ()14p. To remove redundant constraints, a natural 
strategy is to first compute the sequence {(i'(^), then to start from ti and 

choose among all constraints in this subset the one that "goes furthest" (say ^0,5 above). 
Then go to is and do the same; and so on. This allows to skip the redundant constraints. 
Still, this is not quite what we want because there might be nonzero t^'s in the "overlaps" 
of the chosen constraints (in Figure IHJ this would correspond to nonzero t^'s between 
Ti(^h+i) aiid Tj{h))j which again would give strict inequality in ((11]) . It turns out that 
using the above strategy but starting from the end gives the desired properties. The 
result of this will be denoted by {(^(A;),i(A;))}|^]^ below. To see why this works, consider 
Figure El Say that at some point in our choice of saturated constraints we get to toll 
arc rj(/j+2)- We then choose the constraint covering Ti(^h+2) (among {{i' {k) ^ j' {k))} j^^^) 
which reaches furthest in the direction of s, here denoted All arcs between 

Ti{h+2) aiid Ti(^h+i) are now covered (some of those are not represented in Figure EJ, so 
we go to Ti(^h+i)- Likewise, we choose the path Now note that when we set 

the recursion (|lUjl . the constraint corresponding to became saturated 

(or it had been saturated by a previous toll), so that all toll arcs between Tj^/j^^) and 
Tj(^h) (iiot represented here) have been set to 0, as claimed. The heart of the proof that 
follows is to check that this is indeed the case. 
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Consider the following construction. Start with k = 1. In the recursive formula UlUj) . 
let {i' (k) , j' (k)) denote the index for which the minimum is attained when setting tk- 
This makes sense because we assumed that there exists a toll-free path from s to t. In 
case of nonuniqueness, select the largest index j and the corresponding smallest index 
i. Note that there holds 

ti = 0, \f k + 1 <l <j'{k) -1. (15) 
Then rather than evaluating H10|) for tfc+i, we jump from tk to ij'(fc) and set 

ii'{l),j'{l)) = ii'ik),j'ik)), y k + i<l<fik)-i. (16) 

Recursing gives a sequence of indices {{i' (k), j' (k))}/^^^. 
In view of Theorem 0J which implies that 

TiTk) <'Ui^h),j{h) - ^i{h),j{h), (17) 

k=i{h)+l 

we say that covers the toll arcs for z(/i) + 1 < A; < j{h) — 1. To derive ((TUl. 

we look for a subset of {'yj'(A;)j'(fc)}fc=i that covers all toll arcs of P and such that = 
for all arcs Tk covered by more than one subpath. 

We proceed backwards. Select li = and recursively compute Ik = i'{lk~i) until 
/g+i = for some index q. There follows: 

= m,^ + 1 and i{lq) = 0. (18) 

Now, reverse the sequence by setting {i{k), j{k)) = (i'(/q_|_i„fc), j' {Iq+i^k))-, ^ 1^ k < q, 
to obtain a sequence that satisfies the assumptions of Theorem ^ 

1. i(l) = 0, j{q) = + 1: This follows from (UHj). 
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2. i{h) < i{h + 1): This follows from the construction of the backwards sequence 
{k}l=v 

3. i{h + 1) < j{h): By contradiction, assume that j{h) < i{h + 1). This implies 
that rj(ft) is not covered by a toll-free subpath. This is impossible by the very 
construction of {{i{k), j{k))}1^^. 

4. j{h) < i{h + 2): By contradiction, assume that j{h) > i{h + 2). Then 
{i'{lq+i^h), j'i^q+i-h)) = (iih), j{h)) implies that = i{h) for all indices Iq+i^h < 
I < j'i^q+i-h) = j{h) and in particular for index I = Iq+i-h-i = i'{lq+i-h-2) = 
i{h + 2) (< j{h) by assumption). This implies the contradiction = 
i{h + l) = i{h). 

5. (|13() and (|14|) : By construction, every arc is covered by at least one path. By the 
preceding inequalities, the only arcs covered by more than one path must belong 
to the interval k = i{h +!) + !,... ,j{h) — 1, for some index h. By ()15|). their tolls 
must be zero, and is satisfied; H14|) then follows from (jl^f) . 

□ 

3.4 Partitioning the set of toll arcs 

The approximation algorithm progressively removes toll arcs from the network. It makes 
use of the following definition. 

Definition 3 A descendant P' of a valid path P G A/q is a simple path from s to t in 
Nq{P) that traverses a subset of the toll arcs of P in the same order as does P. 



25 



Let P be a valid path in TVq. Theorem [S] suggests a way of constructing a descendant 
of P that stands a chance of achieving a high revenue, whenever P's revenue is low. Let 
us consider the set 

Vi(^2),j{2), • • • , Vi{q),jiq) (19) 

of toll-free paths such that 

= i(l) < i{2) < i(l) < i{3) < i(2) < i(4) < i(3) < . . . < i{q) < j{q-l) < = 

(20) 

and equality 1)14(1 holds. If the maximum revenue on P is B{P) then descendants need 
not be considered, because their upper bounds are smaller than or equal to B{P). This 
happens in particular q = 1. We can therefore assume that q is larger than 1. 

We now consider two descendants: Pi contains all Ui[h),j(h) with h odd and is com- 
posed of arcs of P between them; Pj is constructed in a similar manner, with even 
values of the index h. For instance. Pi may start in s, borrow fj(i) j(i), take path P 
between iNiT(Tj(i)) and TERM(rj(3)), bifurcate on ^^i(3)j(3), return to P, and so on. Such 
a pattern, which is allowed by (120(1 . is performed by procedure TollPartition (see 
Figure 1111)1 . Note that Pi and P2 have no toll arcs in common and that some toll arcs 
of P belong neither to Pi nor to P2. The rationale behind this construction is the rela- 
tionship between the maximum revenue achievable on P and the upper bounds on the 
tolls of Pi and P2 given by Theorem |21 

Both these descendants are valid paths. If this were not the case for Pi, there would 
exist a subpath {h < h' both odd) such that lAn^h-^^ji^y^ is strictly smaller than 

the length of Pi between TERM(Tj(/i)) and iNYT{Tj(^y^) (indeed, a path is valid if and only if 
the null toll vector is consistent with it). Since this length is equal to >Cj(/j)j(/j/)-|-7^(/j)j(/j/) 
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by (|14() we obtain that the toll assignment on P is not consistent, a contradiction. A 



similar argument applies to P2- 


Algorithm TollPartition 




Input: a path P and sequence of indices {{i' (k) , f {k))} j^^^ 




Output: two paths Pi, P2 




• Set I := m^, r := 1 




• (Backwards selection of constraints) While / > 0, do 




- Set (z"(r),/(r)):=(^'(/),/(/)) 




- Set / := r := r + 1 




• Tjet n be the last index in the loon above 




• (Reversing the order of the sequence) For all 1 < < g, set 








• (Constructing the descendants) Let Pi contain all Vi(^h),j{h) with h 


odd and 


be composed of arcs of P between them; let P2 be constructed in 


a similar 


manner, with even values of the index h 




• Return Pi , P2 





Figure 10: Algorithm TollPartition 
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Figure 11: An instance of MaxToll. Arcs without label have cost 0. 
3.5 A detailed example 

We apply the algorithm to the example of Figure [TTl Here the shortest path is of length 
and shortest toll-free path has length 10 so that LP = 10. We start with the shortest 
path when tolls are set to 

Pq = {s, Vi, V2, V3, Vi, V5, V6, Vj, Vs, t). 

We apply MaxRev: Ti = T{ti) = 1 because of {s, V3) (i.e. ^0,2) so (i'(l), j'(l)) = (0, 2); 
T2 =T{t2) = 2 because of (s, v^) (i.e. fo,4) so (i'(2), j'(2)) = (0,4); now we jump over rs, 
setting Tg = r(T3) = and (i'(3), j'(3)) = (0, 4); finally = T{t4) = 1 because of (^4, t) 
(i.e. ^2,5) so (i'(4), j'(4)) = (2,5). The total profit is 4. We then run TollPartition. 
We start from the end, i.e. r4. We get {i" (1) , j" (1)) = (i(4), j(4)) = (2, 5). We jump over 
to T2. We set (i"(2),i"(2)) = (i(2),j(2)) = (0,4). And we are done. We finally obtain 
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the sequence (i(l), = (0,4), (i(2),j(2)) = (2,5) (note in passing that the only toll 
arc in the "overlap" of (f(l), j(l)) and (f(2), j(2)), i.e. rs, has been set to 0, as claimed by 
Theorem inj . The new paths are thus Pi = (s, vy, vg, t) and P2 = {s, vi, V2, ^3, v^, t). 

Now consider Pi. We fix Ti = T2 = T3 = +00. The only constraint on T4 comes 
from the arc (s, t). So we set T4 = 10 — 3 = 7. There are no descendants to this path. 

Consider P2. Fix T3 = T4 = +00. Apply MaxRev again. The toll Ti is bounded 
by 1 because of {s,v-s). Then T2 = 10 — 1 — 1 = 8. The total profit on P2 is 9. Again, 
applying TollPartition, we get that there are no descendants (the descendant is 
actually (s, t)). 

Summing up, the best profit was achived on P2 with 9. Note that in this example, 
the algorithm returns the optimal path, but this is not always the case. 



4 Analysis of the algorithm 
4.1 Approximation bound 

We shall prove that ExploreDescendants is an ^ log2 mx + 1-approximation algo- 
rithm for MaxToll. The exact approximation factor is given by the recursion 

a{k) = - max {1 + a{i) + a{j)} , (21) 

2 i+j<k 
0<i<j<k 

with a(l) = 1. It can be shown by induction that for all k, 

a{k) < ^ log2 k + 1 

Definition 4 The maximum revenue Vp induced by path P is sufficient if 



1 

a(m^) 



Vp > -T-prB{P)- (22) 
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Theorem 6 Let P be a valid path on Mq. If the maximum revenue achievable on P is 
not sufficient, then either path Pi or P2 (say P') returned by TollPartition satisfies 

^ TzB{P') > -^BiP) (23) 



Proof: Let 

with corresponding values of Uij, Ck^i and T^^i for all i, j, k, I, k < I. Similarly, for 
r = 1, 2, 



Pr — (^0,1) ■'"15 '^1,2) ■ ■ ■ ^'^^PrJ '^^Pr ^Pr 



with corresponding values of Ul j, ^ and 7^''^ for all i, j, k, I, k < I. Let 

be the sequence of toll-free paths obtained from Theorem |SI 

For all h > 1, one of the paths Pr contains ^Ji(/i-i).j(/t~i)5 while the other contains 
'^i{h),i{h)- Therefore, the subpath of P between TERM('7:j(;j)) and iNiT(rj(/j_i)) is in neither 

Pi nor P2. The remaining arcs of P belong to either Pi or P2. This implies that 

<? <? 

h=l h=2 



Y [^i{h),j{h) - ^iih),jih)] + 2/:o, 



rriy +1 



h=l 



By Theorem [SI the first term on the right-hand side is equal to Yl^=i and is strictly 
smaller than B{P) / a{mf^) by hypothesis. Multiplying by (—1) and adding '2LIq^p_^^ = 
lA^ p. ^ , on both sides, we get 

0,m,j,-^+l 0,mj,^+l " 

B{Pi)+B{P2)>(2-—^\b{P). 

\ OL [TTlrp J J 
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By definition of a: 

a(m^) > Ml + a("ir ) + «("^T )1 =^ (2 r-pv) 

2 1 J \ a(^myj / 



1 ^ ^ a(my^)+a(my^ 



Q!(?ny) 



By substituting in the preceding inequality, we obtain 

Va(m^') / \a{mp) J \a{mif) J 

which yields the desired result. Indeed, if we had —^B(Pr) < -T^B(P) for r = 1, 2, 
this would imply the opposite inequality. 

□ 

As a corollary, we obtain the main result of this section. 

Corollary 2 Let APP denote the revenue obtained from the application of the procedure 
ExploreDescendants to a path Pq of shortest length Cq in A/q. Then, 

APP > — ^^LP > —^—OPT. (24) 

a{mtj?) a{mT) 

Proof: Let P be a vahd path and (V, T, P) the output of ExploreDescendants (P). 
It is sufficient to show, by induction on that 

V>-^B{P). (25) 

This statement is true if = 1 since the upper bound B{P) is always achievable on 
a path with a single toll arc. Now assume that the property holds when the number of 
toll arcs is less than > 1. Let {Vp, Tp) := MaxRev(P). If Vp is sufficient, ((^ is 
satisfied. If Vp is not sufficient then, by Theorem |21 TollPartition returns a path P' 
with , ^„i, B(P') > , ^ p, B(P). Since the number of toll arcs in P' is less than that 

a(m^ ) a(mf ) ^ ' 

in P, property is satisfied for P', and the preceding inequality implies that is 
satisfied for P as well. 
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□ 



Note that this result apphes to the case of negative tolls as well since the upper 
bound is the same. Indeed, ExploreDescendants allows only nonnegative tolls 
but it computes a feasible solution with a revenue at least LP /a, where LP is unchanged 
in the unbounded case. 

4.2 Tightness of the approximation 

The approximation algorithm determines, in a constructive manner, an upper bound 
a{mT) on the ratio OPT/ LP. In this section, we show through a family of instances 
that this bound is tight. 

Theorem 7 Let T{mx) denote the set of instances of MaxToll correponding to a 
fixed number of toll arcs mx- Then for all niT > 1, the relaxation gap on X{mT) is 
a{mT), that is 



where LP[I] and OPT[I] denote respectively the upper bound and optimal value for 
instance I. 

Proof: Let us consider a two-node and two-arc network Z{1), with origin si and destina- 
tion ^1. The first arc (a toll arc) has fixed cost and the second arc (a toll-free arc) has 
fixed cost 2. We recursively construct a network Z{k) made up of a copy of Z{\k/2'\) and 
a copy of Z{[k/2\) (if k is even, there are two copies of Z{k/2)), and two distinguished 
nodes, the origin Sk and the destination tk- These are linked by five toll-free arcs, as 
illustrated in Figure IT^ Strictly speaking, the vertices S[A:/2] ) ^[k/2'] ) ■^lk/2] ) ^lk/2\ should 




(26) 
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be re-labeled, otherwise many vertices will share the same label even though we consider 
them all distinct. The parameters and bk of Z{k) are set to 



Ofc = 1 + a 



a 



bk = l- a 



+ a 



(27) 



with ai = 61 = 1 and a is defined in (|2T|l . Note that a^, bk > because < a([A;/2]) 
a(LA:/2j)<l. 



It is not difficult to show, by induction, that 



a{k) 



l + a 



+ a 



(28) 



Let LP{k) and OPT{k) denote, respectively, the relaxed and optimal revenue values 
on Z{k). Let hl{k) be the length of the shortest toll-free path from to in Z[k). 
We claim that U{k) = 2a{k). Indeed, we have U{1) = 2a(l) = 2 and, assuming that 
U{k') = 2a{k') for aU k' < k, it follows from ^ that: 



U{k) = min <^ 



+ bk, U 



+ CLk, U 



+ U 



2a{k). (29) 



Clearly, 0PT{1) = LP{1) = 2 and LP{l)/OPT{l) = a(l) = 1. For k > 1, the shortest 
path in Z{k) with tolls set to has length 0. This implies that LP{k) = U{k) = 2a{k). 
To conclude, we need to show (by induction) that OPT{k) = 2. We consider two 
cases. If the optimal path on Z{k) goes through the arc joining i[A:/2j a-^d then 
OPT[k) < 2 because + b^ = 2. We actually know in that case that OPT{k) = 2 
because, by induction, a profit of 2 is achievable inside each of Z{\k/2~\) and Z{\k/2\) 
separately and the only new constraints are those associated to and bk- Otherwise, 
the optimal path belongs entirely to Z{\k/2\) or Z[\k/2~\) in which case OPT{k) = 2 
by the induction hypothesis. ^ 



^It is straightforward to check that the same argument works in the negative tolls case. 
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Figure 12: Networks Z{1), left, and Z{k), right. Toll arcs are dashed. 

□ 

We have just proved the optimality of our approximation factor with respect to the 
upper bound. We can prove more than that: the same family of instances, under a 
slight modification, can be used to show that our analysis of ExploreDescendants 
is tight. For this purpose, we require instances where APP is much smaller than OPT. 
This is not the case in the examples of Figure ^1 where, indeed, APP = OPT = 2. 
However, our aim is attained if we add to Z(k) a toll arc of fixed cost 1 from to t^- 
Then OPT = LP — 1 = 2a{k) — 1 (the optimal path being the new toll arc) yet APP is 
still 2 since the algorithm starts with a path of length zero and misses the optimal path 
of length 1. 

4.3 Running time analysis 

MaxToll is initialized with a shortest path Pq, which can be computed in O(n^) time. 
The toll arcs of the descendants constitute a subset of the toll arcs in Pq, and their 
traversal order is the same. Therefore, the values Uij, i < j computed for Pq can 
be reused for all descendants, under an appropriate renumbering. This operation is 
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achieved in 0{'m^^n?) = 0{m-Y'n?) time. The time required to evaluate C-ij, i < j, as 
the algorithm proceeds, is much smaller. 

Within ExploreDescendants, MaxRev requires at most 0{{m^)^) to compute 
the maximal revenue induced by path P. Based on the indices {(i'(^), ob- 
tained from MaxRev, TollPartition generates two descendants in 0{m^) time. It 
follows that the running time of ExploreDescendants on a path P is determined by 
the recursion 

T(m^) = T(m^i) + T{m^^) + 0{{m^f) (30) 

where + < m^. Therefore, the worst-case complexity, achieved when 

is always equal to — 1, is 0((m|?)^). The worst-case running-time of the entire 

algorithm is 0{mT{iTij. + n'^)). 

5 Concluding remarks 

Our algorithm can also be applied to the multi-commodity extension of MaxToll con- 
sidered in ^ni, where each commodity k £ IC is associated with an origin-destination 
pair. Given a demand matrix, users solve shortest path problems parameterized by the 
toll vector T. If distinct tolls T'^ could be assigned to distinct commodities, the multi- 
commodity extension would reduce to a |/C|-fold version of the basic problem. Otherwise, 
the interaction between commodity flows on the arcs of a common transportation net- 
work complicates the problem, both from a theoretical and algorithmical point of view. 
Of course, we can obtain an 0{\JC\logmT) guarantee by applying MaxToll to each 
commodity separately and then selecting, among the |/C| commodity toll vectors, the 



35 



one that generates the highest revenue. However bad this bound is, we conjecture that 
it is tight with respect to the relaxation, which is the sum of the single-commodity 
bounds, weighted by their respective demands. Indeed, we believe that the instances of 
Figure IT^ can be generalized to the multi-commodity case. 

Other generalizations of MaxToll involve capacity constraints and lower bounds 
on tolls. In the latter case, the relaxation gap becomes infinite for any value of m^, 
and our approach fails, as procedure ExploreDescendants becomes irrelevant. A 
completely different line of attack is then required. 

Finally, we raise the following important issue: Can our ^ log2 rriT + 1 guarantee be 
improved? Such results would obviously require a tighter upper bound than the one 
used in this paper. 
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